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Abstract 

The evolution problem for a quantum particle confined in a f D box and interacting with 
one fixed point through a time dependent point interaction is considered. Under suitable 
assumptions of regularity for the time profile of the Hamiltonian, we prove the existence of 
strict solutions to the corresponding Schrodinger equation. The result is used to discuss the 
stability and the steady-state local controllability of the wavefunction when the strenght of 
the interaction is used as a control parameter. 

1 Introduction 

We consider the evolution problem for a quantum particle confined in a ID bounded domain and 
interacting with one fixed point through a delta shaped potential whose strength varies in time. 
The Hamiltonian associated to this class of potentials, denoted in the following with i?Q,((), is 
defined in terms of a time dependent parameter, a(i), which describes the time profile of the 
interaction. Under suitable regularity assumptions, we study the dynamical system defined by 
Ha[t), its stability properties and the local controllability of the dynamics when a is used as a 
control function. 

General conditions for the solution to the quantum evolution equation related to non au- 
tonomous Hamiltonians, H{t), have been long time investigated (e.g. in [5D], [15] and [H]). 
When the operator's domain D{H{t)) depends on time (as in the case of time dependent point 
interactions), the Cauchy problem 

" = -^H{t)i, 

was explicitely considered in [15] by Kisynski using coercivity and Cf^^^-regularity of t ^ H{t). 
A similar assumption, a € Cfg^{tQ ^-\-oo) in the above notation, is used in the works of Yafaev 
[3T], [35] and [53] (with M. Sayapova) to prove the existence of a strongly differentiable time 
propagator for the scattering problem with time dependent delta interactions in R"^. Such a 
condition, however, can be relaxed by exploiting the explicit character of point interactions. The 
dynamics associated to this class of operators, indeed, is essentially described by the evolution 
of a finite dimensional variable related to the values taken by the regular part of the state in 
the interaction points (e.g. in [TT], [5]). This leads to simplified evolution equations allowing 
rather explicit estimates under Fourier transform. Using this approach, the quantum evolution 
problem for a ID time dependent delta interaction has been considered in [14] . In particular, it 
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is shown that a G Hi{0,T) ahows to define a strongly continuous dynamical system in iy^(M), 
while a G iJ4(0,T) leads to the existence of strong solutions to 



{ 



ipt=o = V'o e D{Ha{0)) 



(2) 



The same strategy has also been adopted to study the difi^usion problem with nonautonomous 
delta potentials in M.^ [11) . and the quantum evolution problem for a ID nonlinear model where 
the parameter a is assigned as a function of the particle's state [21 • 

The techniques used in these works can be adapted to the ID confining case. However, the 
lack of a simple explicit expression for the free propagator kernel and the use of eigenfunction 
expansions, which replaces the Fourier transform analysis in this setting, necessarily requires some 
additional efforts to obtain a result. In the perspective of the stability and controllability analysis, 
we restrict ourselfs to the case a G H^{0,T) and: Ha = + cw5 with Dirichlet conditions 

on the boundary of / = [— 7r,7r]. Under these assumptions, we prove that the evolution problem 
^ admits a strongly differentiable time propagator preserving the regularity of the initial state. 
This turns out to be a key point to study the controllability of the dynamics. 

The controllability of a quantym dynamics through an external field has attracted an in- 
creasing attention due to possible applications in nuclear magnetic resonance, laser spectroscopy, 
photochemestry and quantum information. This problem has been considered for confining 
Schrodinger operators with regular potentials of the form: H = — + V{x) + u{t)W{x), where 
is the Dirichlet Laplacian in the bounded domain Vt C M™, while u is used as a control 
function. The particular setting: m = 1, V — Q and W{x) — x, corresponding to a quantum 
particle confined in a ID box and moving under the action of a time dependent uniform eletric 
field, has been considered in [5]. The exact controllability of the quantum state was proved, 
in this case, in i?^-neighbourhoods of the steady states by using u G L^{Q,T) controls with 
T > Tm > (actually a simpler version of this proof holding for all T > is given in 8^ ) . For 
the same system, the exact controllability between neighbourhoods of any couple of eigenstates 
is discussed in 7 . In the more general setting, with V,W £ C°° (f2,R) and m being any space 
dimension, differents approximate controllability results in have been presented in [101 and 
|16) . An example concerned with singular pointwise potentials is presented in jT|, where point 
interaction Hamiltonians are used to construct a general scheme allowing to steer the system 
between the eigenstates of a drift operator Hq (the ID Dirichlet Laplacian). The idea is to in- 
troduce adiabatic perturbations of the spectrum of Hq producing eigenvalues intersection, while 
controlling the state's evolution through the adiabatic theorem. The result is the approximate 
state to state controllability in infinite time. The particular case of ID point interactions is also 
related to control problems on quantum graphs where these Hamiltonians naturally arise. 

In Section 31 the nonlinear map a — )> Vt, associating to the function a the solution at 
time T of is considered for a fixed initial state ipo G £'(iJct(o))- This can be regarded as a 
control system where the control parameter is a and the target state is ipT- Using the regularity 
properties of the time propagator, we show that: for -00 € H Hq{I), the map a — ?► i/jt 
is of class C^(-ffQ(0,r), n Hq{I)). Then, the local controllability relies on the surjectivity 
of the corresponding linearized map according to an inverse function argument. This point is 
considered in Section 14.21 where a general condition for the solution of the linearized control 
problem is given in (|98)) - (|99)) and proposition [161 When ipQ coincides with an even eigenstate of 
the Diriclet Laplacian, this scheme provides with a result of local steady state controllability in 
finite time. 
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2 The Model 



Point interactions form a particular class of singular perturbations of the Laplacian supported 
by finite set of points. In W^, d < 3, these Hamiltonians have been rigorously defined using 
the theory of selfadjoint extensions of symmetric operators (e.g. in @]). The definition easily 
extends to the case of bounded regions by taking into account the Dirichlet conditions on the 
boundary for the functions in the operator's domain (e.g. in [5]). In what follows we consider a 
delta perturbation of the Dirichlet Laplacian, centered in the origin of the interval / = [— 7r,7r]. 
In terms of quadratic forms, this family of Hamiltonians acts on Hq{I) as 

H^ = -Af+aS (3) 

where Af is the Dirichlet Laplacian on the interval, S is the Dirac distribution centered in 
the origin, while a is a real parameter. The operator's domain D{Ha) extends to all those 
vectors ip S L^{I), such that: Haip G L'^{I)- The description of D{Ha) is strictly related to the 
properties of the Green's kernel, Qq, for the resolvent of the unperturbed operator (— + z) 
Fix x' e /; whenever z does not belong to the spectrum of A^^, Qq satisfies the equation 

{-Af + z)g^{x,x')^6ix-x'), zeC\(j^n (4) 

The solution to ^ can be represented as the sum of the 'free' Green's function plus an additional 
term taking into account the conditions at the boundary 

^-■/^\x-x'\ 

g^{x,x') = ^-^ h{x,x\z) (5) 



{-Af + z)h{-,x',z) = 



(6) 

X — ztlT 



The solution of ^ is 



(7) 

Another useful representation of C/q is in terms of Fourier expansion w.r.t. the Laplacian's 
eigenfunctions. Denoting with ipk and respectively the eigenfunctions and the eigenvalues of 
—Af, we have 

{-j= sin ^x k even 
^cosf. kodd ' ^^ = r '^"^ 

and 

eo^(x,a;') = ^-^Vfe(a;')V'fc(a:), z^<yAf (9) 

Using ([3]) and (|3]), a straightforward calculation shows that Hatp G L'^il) whenever ip has the 
form 

V'-0 + 9^o(-,O), 0ei?(Af), -g^aV^O) (10) 

This is a general characterization of D{Ha). It can be derived by using the von Neumann theory 
of selfadjoint extensions to identify Ha as the extension of the symmetric operator Hq 

D{Ho) = yjeH^nH^{i)\m = o} 
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to those vectors £ D{H^) fulfilling the selfadjoint boundary conditions (e.g. in [5] and |18j ) 

■ijj{0+) -TpiO-) = 



i^'(0+)-V'(0-) = ai^(0) ^^^^ 

The following proposition is a rephrasing of the result in in the bounded domain case. 

Proposition 1 Let a G M. A G C\(t^d and denote with Ha the family of selfadjoint extensions 
of Hq associated to the boundary condition The following representation holds: 

D{Ha) = G {0}) n H^il) 1^-0^ + qQ^ (-,0) ; (f,^ e n H^iiy, -q = a^{0) } 

(13) 

Haij = ^^4>^-\qg^{;0) (14) 
Moreover, for z G C\R, ip G I? (J), the resolvent {Ha + z) ^ writes as: 

{He. + zy' ^ = (-Af + zy\^ i + ag^{o) ("^^ + "^^^^ ^° ^^^^ 

Remark 2 Although the representation: = (f)^ + qQ^ (-,0) ma?/ change with different choices 
of X, the operator Ha depends only on the value of the parameter a related to the interaction's 
strength. 



Ha(t) 



Next we assign a as a function of time and consider the non autonomous system defined by 

(ifip{x,t)=Ha(t)i'{x,t), 

, at ^^g^ 
The mild solutions are 

V,(., t) = e**^? i^o + t f q{s) e''^'-'''^" S ds , (17) 



where e**'^? is the time propagator associated to — Af and —q{t) = a(t)V'(0, t) fixes the boundary 
condition of the operator's domain. The action of the operator e**^^ on L'^{I) is 

e^*^"/ = ^(V'fc,/)i2(,) e-^^^*Vfc yfeL^I) (18) 

fcSN 

/TT 
uv. This relation 
-TT 

suggests to replace the distributional part at the r.h.s. of (fT7)) with 



^ V f q{s)e-'^-'^'-'Usi^k 



It follows ^ 

V'(-,t) =e"^?^o + ^ V / <?(s)e^*^'=(*-^)dsVfc (19) 



fc odd 
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q{t) = -a{t) 



7^ tTi. Jo 



k odd 



(20) 



The essential informations about the dynamics described in p9p - p0p are contained in the aux- 
ihary variable q, usually referred to as the charge of the particle (e.g. in [3]). In what follows, 
the solution of the above problem is considered under suitable regularity assumptions for a. The 
equivalence with the original Cauchy problem ([TB)) is further addressed. Finally, the dependence 
of this solution from a, considered as a control parameter, is investigated. Our result is exposed 
in the following theorem. 

Theorem 3 Assume: a £ H^(0,T), ipo £ D{Ha(o)) and let V'T(^Oja), q{t, ipo, a) respectively 
denote the solution at time T and the charge as functions of the intial state and of the parameter 
a. The following properties hold. 

1) The system (W^ - admits an unique solution ipt S C{0,T; ffi(/)) n (0, T; L'^{I)) with: 
ipt & ^(-f^Q(t)) o-iT-d idtipi'it) — H^(^t-fTp{-,t) at each t. 

2) The map a ■i/;T(^OjQ^) *s C'^(i/^(0, T), iJg (/)) in the sense of Frechet. Moreover, the 
'regular part' of the solution at time T, defined by: tpT{'4'0i<^) ~ Q(T,'4'o,ct)QQ and considered as 
a function of a, is of class C^{H^{Q,T), H"^ r\HQ{I)). //a G -ffg (0, T), then •iljT{'4'Q,Oi) coincides 
with its regular part and a Vt(V'o,Q;) G C^{Hq{0,T), H"^ n Hl{I)). 

3) Let a G _ffQ(0,r) and W he the subspace of nHQ{I) generated by the system: {V'fc, k odd}. 
If ipo e W, then: a Vt('0o, a) e C^{H^{0, T), W). // Vo = ipk for a fixed k odd and T > 8n, 
then: it exist an open neighbourhood V x P C Hq{0,T) x W of the point (0,?/'^) such that 

, Q!)|y : — V P is surjective. 

The proof is developed in the Propositions [T^ [HI and in the Theorem [T71 of sections 3 and 

4. 

Notation: The Sobolev spaces are denoted with H'^ or Hq in the case of Dirichlet boundary 
conditions. The notation '<' is to be intended as '< C" where C is a suitable positive constant. 



3 The existence of the dynamics 

In what follows the existence of solutions to (|16p is discussed under the assumption: a G (0, T). 
Our strategy consists in using the mild approach described by (|T5)) - (P(I)) . It articulates in two 
steps: First we prove that (j20l) admits an unique solution in H^{0,T), for any finite time T and 
any a G H^{0,T). This result is then used to obtain estimates for the solution of (fT6|) . 

3.1 The charge equation 

Let consider the linear map U 

ft 1,2 

Uq=^ q{s)e-'^>'^'-''>ds, D^{keN, k odd} , Afc = — (21) 

Remark 4 Depending on the value of T , some of the functions e"'^''' may belongs to the stan- 
dard basis {e'""*; w = ^j^^^ of the space L'^{0,T). In particular, for: T = SirN, N e N, it 
follows 

K''^'}.eMc|e'-"*;u^ = ^| (22) 
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The i22\) will be often used as an auxiliary condition in the forthcoming analysis. 
Lemma 5 Let Ht denotes the Hilbert subspace 

HT^{qeH\0,T), qiO)^0} (23) 
equipped with the H^-norm. The map U i21\) is bounded in Ht with the estimate 

\\Uq\\n^<(T+\)\\q\\^^ (24) 

Proof. We use the standard basis |e'""*; u = ^}„g2 -^^(0'^)- The Fourier coefBcients of 
g(s)e^''*''=*^*~'*-'ds w.r.t. the vectors e*""* write as 

dt [ dsq{s)e~'^^^''''>e''^"* = [ dsq{s)e'^^'[ dt e-'^^^+^'''>\ 

Jo Jo J s 

Assume that T fulfills the condition ([^ and, consequently, {e"*'*'''*} forms a subset of {e"^"*; w = 
This choice allows to avoid small divisors, while it does not imply any loss of generality. With 
this assumption, the above expression writes as 



I dt / dsq{s)e 
Jo 



/J" ds q(s)e'^>'' {T - s) , for At + w?i = , 

/g^ ds q{s)e'^>''' - ds q(s)e-*'^"*] , for Xk + ujn ^ . 



Xk+ujn 



and the expansion 



T pt 

dt / dsg(s)e-*^'=(*-'')e-*'^"* = 
Jo 

= {Tq^Nk^^itq)^r^,.)e-^^^'+ V ^^[g_^,.-g„]e-"* (25) 

Xk+i^n^O 

follows, qn denoting the n-th Fourier coefficient of g (in particular: q-Nk^ = (^i '2~*^''*)l2(o t))' 
Let S* = {n e Z : wn ^ {^^k}keD }' have 

~ iT 
(te-'^"*) e~"^"* ^ , VneS. 



This allows to write 

Afe + ujn \ I 



where lis is the projector over the subspace spanned by {e*'^"*}„gg- Using ([25]) and the above 
relation, the Fourier expansion of Uq writes as 



Uq^Y. (T,-Nk^ ~ itqL,,.) e^^^^* + E 1-->^^^S -EE at!^ ■ 



(26) 
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Consider the contributions to Uq: 

i) The norm of the first term at the r.h.s. of (PB)) is bounded by 



keD 



<2T|lq|l^. 



(0,T) 



(27) 



L2(0,T) 



ii) For the second term, one has 



keD 



keD 



L2(0,T) 
— £ 

T 



iXk t 



- \l-Nk^f < 1191^2(0^7) 

L2(0,T) keD 



(28) 



ill) The remaining term is a superposition of L^-functions: X^ngs ^fc+Lw e*""* parametrized by 
the index k E D. Its n-th Fourier coefficient is formahy expressed by: X^fcgD xl'+un ''"^'^ 
L^-norm is 



EE 



keD nes 

To study the sum at the r.h.s, the relation 



-E 

L2(0,T) 



E 



-Et^ = ^o(0'0) 



arising from ([5]) and dH]), is used. According to ([7]) and ([50)1 . it follows 
111 ^ 



(29) 



(30) 



E 



This leads to 




r2 sinh (TTVwn) , 

for n > 0, n e S* 
for n < 0, n £ S 



neS. (31) 



(32) 



Consider the map T : ^ fn (j2keD At+c^w ) ^^{S). Using (15^ . T can be identified with 
the limit of the sequence of finite rank maps: T^v/n ^ \ fn { 'l2keD \ +un ) i > ^2(5*)- 

L V ^ J J n<N 

operator norm. Thus, it defines a compact operator in ^2(5') (e.g. in Theorem VI. 12, JJj), and 

we get 



E 

nes 



E 

keD 



Xk + ujn 



E 

\keD 



qn 



Xk + ojn 



\Tiqn)\\ 



< 



l^2(S) 



£2(S) 



2^ X,. 4- / 



keD nes 



^qn ^iumt 

Xk + u!n 



< 



(33) 



(34) 



L2(o,T) 
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The estimates ^ and (|M1) yield 

\\Uqh 2(0,T) S{T +1) lkllL2(0,T) 



(35) 



This can be obviously extended to any finite time T E (0, +oo) by fixing Nrp G N such that: 
SirNf > T and setting T = SttNj, in ((2^ . 

Next, consider q e Ht- An integration by part of (1^ gives 



t^D 



(36) 



Due to the asymptotic behaviour of the coefficients ^ p- , this sum uniformly converges to a 
continuous function of t € [0, T], with 



From ([55]) . we get 



and 



< 



(7^ + 1)11911^2(0,7) 



L2(o,T) 



(37) 

(38) 
(39) 



In the next Lemma we give an iteration scheme for the solution of charge-like equations. 

Lemma 6 Let f,ip E H^{0,T), ip real valued, and Gq{',0) be the Green' function defined in 
with X E C\R. The equation 



v = f-ip{v{0) e''^' g^{;0) +-Uv 

\ x=0 TT 

has a unique solution in H^{0,T) allowing the estimate 

II"IIhi(o,t) ^ II/IIhi(o,t) (l + II'<^IIhi(o,t) {\\v\\h^q,t) 
where V{-) is a suitable positive polynomial. 



(40) 



(41) 



Proof. Since Lemma [5] provides with an estimate of U f in the "HT-norm, we introduce the 
variable w{t) ~ v(t) — ^{0). In this setting, our problem writes as 



w — R — ip—Uw , 



Rit) - fit) - viO) - ipit) v{0) e^'^' ^o^(-,0) 



+ -u[vm (t) 

x=0 TT 



(42) 



(43) 



Next, the regularity of R{t) is considered. Using the definition of U, and the explicit expansion 



TT ' * 



keD 



Afc + A 



(44) 



R{t) writes as 

Rit) = fit) - vio) - vioMt) l^e**^? g^i;0) 

= fit)~viO)~viOMt)-Y,(^-e 



keD 



-iXkt 



(45) 



where the last term at the r.h.s. is H^{M.) for any —A ^ f-Af ■ RecaUing that H^{0,T) is an 
algebra w.r.t. the product, the above relation and f,ipG H^{0,T) imply: R S H^{0,T). Taking 
into account the initial condition R{0) — (following from the definition of the rescaled variable 
and the operator U in (gl])), we get: Re Ht (see (^5]) ). 

Next consider the term ip^Uw in (H^ . For w E Ht and (p e H^{0,T), the Lemma [S] yields 



(fi—Uw 



■Ht 



<C^-^MHHo,t)Mnt 



(46) 



for any t <E [0, T] and for a suitable C > 0. To construct a global solution, we notice that, for 



fixed T,C > 0, if £ H^{0,T), there exist S > and a finite set {ij}^]^, io = and 



fulfilling the conditions 



Cil + S) \MHi(t,^,.t,) <'^C5 foraUj, 



(47) 
(48) 



with C5 < 1, depending on 5. According to (^51) and pS)) with j = 1, fyS^C/ is contraction map in 
"Ht^ and the equation (j42p - (l43p admits an unique solution, wi = f [q f^jw, bounded by 

Using dini), the initial value of v if formally related to /(O) by: u(0) (l + ^(0)^^(0,0)) = /(O); 
this allows to define v{0) privided that (f + (p{0)Qq{0, 0)) is not null. According to the resolvent's 
formula ([T5|) . the condition: (1 + ip{0)Qq{0, 0)) — 0, with z ^ C-Af j corresponds to the eigenvalue 
equation for the Hamiltonian H^. Since this is a selfadjoint model, the assumption A S C\M 

implies: (1 + ^(0)^^(0,0)) ^ 0. It follows that: u(0) = (l + (^5(0)^^(0, 0))" V(0) and |w(0)| < 
1/(0)1 .This leads to 

W^lW-Ht, ^ li/llHi(0,fi) (l + ll'^llffMO.ti)) ■ 

The definition: wi = l[o,ti] (^ — ''^(0)), provides with a similar estimate for the function v 



This can be extended to larger times by the following iteration scheme: 

Wj+i{t) = v{t)-v{tj), te[tj,tj+i], j = i,...,iv-i 



(49) 



(50) 



Wj+iit' + t,)^R,+iit' +t,)-ipit' +t,)-[Uw,+ii-+t,)]it'), i'e [0,tj+i-tj] (51) 

TT 
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= fit) - v{t,) ^{t) Uq) e^''^?g^{;Q) 

\ X- 

-y{t) (u [v{t,)] + ^ p (vis) - vit,)) e-^^-^'-^Us] (52) 



I 

— i 



keD 



with the source term Rj+i depending at each step on the past solution l[o,t.](7. From this 
definition, it follows: Rj+i{tj) ~ 0, while an integration by part gives 

^ g (t^ . .,0, (^) - J- f 

(53) 

where have been taken into account. If u e H^{0,tj), all the contributions at the r.h.s. 
of (1551) are H^{0,T) and Rj+i G T-Ltj+i-tj- Then, the contractivity property of the operator, 
expressed by (|5ip . implies the existence of an unique solution Wj+i € 'Htj^i-tj with the bound 

Starting from <\4Q\ . an induction argument based on the itarated use of ((54|) leads to 

II"IIhmo,t) ^ II/IIhmo,t) (i + II^IIhmo,t) +^ (|I<^IIhmo,t))) (55) 

where ■p(-) is a positive polynomial. ■ 

A first application of this Lemma gives an iJ^-bound for the solutions to the charge equation. 

Corollary 7 Let a £ -ff^(0,T;R) and ipo e D[Ha[Q)) and denote with (j)^ the regular part of ipo 
defined according to the representation I113\) . Then: 



i) The equation i20\) admits an unique solution q G H^{0,T) such that 

Mh^o^t) <C'A,0i^l|a|lHi(o,T)(l+'^(ll"llffi(o,T))) (56) 

where V (•) is a positive polynomial, is a positive constant depending on X E C\M and (f)^. 

ii) For a fixed "00 G D(Ha(Q)), the map a q is locally Lipschitzian in H^(0,T). 
Proof, i) With the notation introduced in (I^TI) . the charge equation is 

g==-ae'*'^?Vo(0)-ck-J7g. (57) 

TT 

Using the decomposition: i/'o = 0o + 9(0)^/0 (': 0)> (po ^ <^ H^{I), holding for A e C\R (we 
refer to Proposition [1]), this equation rephrases as follows 

q = _«e**A? _ ^ (^q(0)e«A? g^{.,0)l^^ + ^Uq^ . (58) 

Consider e'*^^ (j)^{0). As a consequence of the Stone's Theorem (e.g. in |2Q])), the operator e* ^ 
defines a continuous flow on n Hq{I) strongly differentiable in L^(/). For 0q e n Hl{I) , 
we use the Fourier expansion (t>^ = ^ aki^k, where the coefficients Ok — {(p^^'tpk) ]^2ii\ ^re 

characterized by: ak G ^2(N), \k0.k G ^2(N) (as it follows integrating by parts and exploiting the 
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boundary conditions of ^fc). The time propagator e'*^? acts like: e'*^?0Q = ^ a^e In 



particular, we have 



fceN 



(59) 



keD 



According to: Xkdk € ^2, it follows: 9(6* ^ 0o(O) = ^*Sfee_D i^fcAfce"*^''*, in the weak sense and: 
gitAf ^A^-Q-j g H^{0,T). Then the first statement follows as an application of Lemma H] with 
/ = -ae"'^?(?!)o(0) and (p = a. 

a) Let a, a € H^{0,T), ipo and ^/iq defined with the same regular part 

^o^^^ + qiO)g^{;0), ,Po^(f>^ + mGoi;0). (60) 

and consider the corresponding solutions to ([57)1 q and q. The initial values depends on a and a 
through the conditions 



9(0) = 



0o'(O) 



1 + ^(0)^0^(0,0) ' 
while the difference u = q — q solves the equation 



m = 



4>m 



1 + ^(0)^0^(0,0) 



S{t) = -{a-a) 



e ^ 00 (0) - - 
a 



(61) 



(62) 



(63) 



From the previous point, we have: e**'^^^ (/)o (0), q, q € H^{0, T). The same inclusion holds for the 
ratios ^, according to the equation's structure. An application of Lemma |5] with: f = S and 
= a, gives 

(64) 



'"llifi(0,T) ^ 11"^ '^IIhi(0,T) 



3.2 Solution of the evolution problem 

Next we consider the system (fT9|) - ([20|) with 

aei/i(0,r), 7Aoei?(ifa(o))- (65) 

Due to the result of Lemma IH the equation (l20l) has a unique solution q G H^{Q,T). This 
can be used to show that the corresponding evolution, defined by p^ . is C([0,r] ,Hl{l)) n 
C^([0,T] ,L'^{I)) and solves the Cauchy problem ([TC)) . Fix A G C\(t^d; from the definition of 
the operator's domain (see P^ *). any i/; e D{Ha(t)) can be represented as the sum of a 'regular' 
part plus a Green's function 

ij{;t)^^\;t)+q{t)g^ (66) 
with q fulfilling (l20l) . In the mild sense, the evolution of the regular part is 

0\x,O=e'*^?^o(a:) + ^ V f\is)e~'^^^'-^UsM^)-qit)Goi^)- (67) 

Next we consider the operator 

F{q, t)^^Y.{f q{s)e-'^'^'-''> ds) ■ (68) 
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Lemmas For q e H^{0,T), the map F{q,t) is bounded m C {[Q,T] , Hl{I)) . 

Proof. Let qt — l(o,t)'? and Nt be the smallest integer such that: SttNt > T. By definition (see 
(mi), {e~*^'=^}j,gpj forms a subset of the standard basis in L'^{Q,%'kNt)- The Fourier coefficients 
of qt along these frequencies are 



/ qt{s)e'^''' ds = / q{s)e'^-' ds = Ck{t) 
Jo Jo 

|cfe(t)|^ < \\qt\\l2(o^s7rNT) = / ds 

ken ° 



and the inequality 



holds. It follows 



and 



sup \\F{q,t)\\l2m < - / k(s)l^ ds^- ||<7||i2(o.r) • 
For q G H^{0,T), an integration by part gives 

F{q, t)^^Y.T- f^W - 9(0)6-*^^* - r g'(s)e-^^'^(*-^) ds) i^k , 
Afe V Jo ) 



(69) 
(70) 
(71) 
(72) 

(73) 



while, using the relation: -r-^'fc = ^k''Pk, the weak derivative 4^F{q,t) can be written as 



d 
dx 



F{q, t) = ^Y.^ f^W - liO)^-'^"' - f 9'(s)e-^^'=(*-^) ds] ■ 



(74) 



For the first term at the r.h.s. of ([71)1 . the asymptotic behaviour of the coefficients A| = | and 

the assumptiom q G H\0,T) implies: E^gd^^ (g(t) - g(0)e-'^'=*) ^/-fc e C(([0, T] , L2(/))) . 
For the remaining term, we proceed as before: let q'^ — l(o,t)9' and denote with c'j. the Fourier 
coefficients of q[ along the freqencies 



4 = / q'(s)e-*^'=(*"")ds. 



It follows 



and 



ken Af. 



sup 

iG[0,T] 



iTT, A? ^0 



L2(/) fce-D 

2 



trk Jo 



-i\k{t-s) 



ds Ipk 



keD 



< h'W 



L2(0,T) • 



LHI) 



The above estimates lead us to 



sup \\F{q,t)\\^,f^i) < MhHo.T) 

te[o,T] ' ' ' 



(75) 



(76) 



(77) 
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Proposition 9 Let a G H^{0,T) and 4>o G ^(-ffa(o))- The system U9\) - i2U\) . has an unique 
solution tpt<^C ([0,r] , 7Ji(/))nCi([0,r] , L^{I)) such that: G D{H^^t)) andidt^pt = H^(t)i^t 
at each t. 

Proof. The proof articulates in two steps. We first consider the conditions i/jj G C ( [0, T] , iJg (/)) 
and ipt e C^([0,T] , L'^{I)). Then we discuss the equivalence of the system (ITll)-([ini) with the 
initial problem. ■ 

Proof. 1) Using the notation introduced in ([68|. the solution ipt of (fT9|) - ((20)l writes as 

V't = e'*^?Vo + i^(g,0 (78) 

with q e H^{0,T) solving the charge equation (pn|) . and V'o G D{Ha(o))- Due to the domain's 
structure, ^po & Hq{I): in this case, the term e^'^^^Vo defines a C([0,T], Hq{I)) map (as a 
consequence of the Stone's theorem, e.g. in [lO])- Moreover, following the result of Lemma |S1 
one has: F{q,t) e C([0,T], H^{I)). In order to study the Ci([0, T] , regularity of V-t, one 

can use the decomposition stated in Proposition[TJ ^po = 4>o+q{0)Go{-, 0), with 0o ^ 
and A G C\cr^p. For w{t) — q{t) ~ g(0), ipt writes as 

i^t^e''^PcP^+F{w,t) + Z{t), (79) 

z{t) = q{o) [e^'^?g>;{;0) + F{i,t) . 

Due to the regularity of (pQ, one has: e'*'^?(/>o e C^([0,T] , L^(/)). Next consider F{w,t): an 
integration by part gives 

F(^,t) = ^ V ^ fw(t) - r w'(s)e-'^''('-^) rfs) Vfe 

and ^ 

^F{w,t) = ^y^[ w'{s)e-'^-^'-'Us4'k^F{w',t). (80) 

Using the estimate ([7^ in LemmalU this relation gives: F{w, t) e C^([0, T] , L'^{I)). Concerning 
the term Z{t), the formula 



ipk 

ken + ^ 
and an explicit computation lead to 

^/ X 9(0) 1 , g(0) Ae-*^"* , 

It follows: Z{t) e Ci([0,r], 7Ji(/)) and: Vt e C([0,r], iJi(/)) n Ci([0, T] , L2(/)). 
2) The regular part of ipt writes as 

c^^ = e"^?^o + F{q, t) - q{t)g^{; 0) (82) 

Using once more the decomposition: ^pQ — cp^ + q{0)QQ{-, 0), cpQ G H"^ n Hq{I), A € C\cr^p and 
w{t) = q{t) - g(0), we get 

0^ = e"^?0o^ +F(u;,i) -«;(i)ao'(-,0) +q(0)g(i) (83) 
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QW= [(e''*^?-l)c?o\-,0) + F(l,i)] (84) 



The regularity of (jj^, yields: e**'^?(/)o e C([0,T] , n -ffo(/)), while, according to the explicit 
form of e'*^?g^{-,0) and F{l,t), Q{t) is 

This yields g(i) G C([0,r] , TJ^ p For the remaining term, F{w,t) - w{t)g^{-,()), it 

has already been noticed that F{'w,t) is continuously embedded into Hq{I); the same holds for 
w(t)5g (-,0), since Qq (-,0) G Hq{I) and w is continuous. Moreover, a direct computation of the 
Fourier coefficients of ^ [F{q,t) — q{t)QQ (-,0)] w.r.t. the system {V'fcjfegN gives 

6, = I ;fe/o^'We-*'^^*"^^rf^ - A«;(t) (t^oM', 0) , ^,)^,(,) fc odd ^^^^ 
1 otherwise 

from which it follows 

- ^ [^^(u;, - u;(i)eo^ (-,0)] - iF{w', t) + \w{t)g^, (•, 0) . (87) 

Then, the estimate ^ in Lemma [8] yield: ^ [F(g, t) - (j(t)^;f} (•, 0)] G C([0,T], L2{I)) and 
[i^(i(;, t) - w{t)g^ (-,0)] e C ([0, T], H^n H^il)) , which implies 



ec([0,r],ij2ni7i(/)) 



this, together with the boundary condition —q{t) = a{t)ijjt{0) (arising from the equation (PU]) ). 
assures that tpt & D{Ha[t)) at each t. 

Next we discuss the last point of the Proposition. The continuity in time of the map ipt and 
(fT9)) allow to write: i/'('ii = 0) = -00- According to ([79| - ((8T|) . the derivative idtipt is 



^J^^t = -^e^'''m+^F{w',t)+^Z'{t) 
, q(0) Ae-'^--* , 



On the other hand, from definition one has 

d2 



J^awV-t = -^<^^ - ^'?(*)^o (-,0) . (89) 



Making use of (|H31)-(Eni) and this writes as 
with 



Hait)A = -T^e'*^^ + tF{w\ t) + Aq(0)go^ (•, 0) - q{{)) — Q{t) 



-liO)-^Qit) = ^FW, t) + ^9(0)^0^ (•, 0) . 



Combining the above relations, it follows 

d 



i-4't - H^(t)H^,t) ^ in L'{I) 
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4 Stability and local controllability 



The evolution problem (jl6p defines a map T(a, ipo) associating to the coupling parameter, a, and 
the initial state, ipo, the solution at time T. According to the notation introduced in (I2ip and 
it writes as 

r{a)^e'^^^o + F{V{a),T) , (90) 
y(a) = g : g(t) = -a{t) (^e^'"^? M^) + '-Uq{t) ) • (91) 

The local controllability in time T of the dynamics described by (jl6l) . when a is considered as 
a control parameter, is connected with the following question: for ipQ € D^H^ij^-^), does a exists 
such that r(a,'0o) = V't for arbitrary V't in a neighbourhood of ^o? In what follows, we fix the 
initial state ?Ao and discuss the regularity of the map a — >■ r(a, V'o)- Hence, this control problem 
will be considered in the local setting where ■00 concides with a steady state of the Dirichlet 
Laplacian and a is small in -ffp (0, T). To simplify the notation, we use r(a) instead of r(a, il^a). 



4.1 Regularity of F 

Recall a standard result in the theory of differential calculus in Banach spaces. Let X and Y 
denote two Banach spaces, W an open subset of X , daF and d'^F respectively the Frechet and 
the Gateaux derivatives oi F : W ^ Y evaluated in the point a. A functional F : F is of 

class if the map 

F' -.W ^ C{X, Y), F'{a) = d„F (92) 

is continuous. 

Theorem 10 ((Theorem 1.9 in [3])) Suppose F -.W ^ Y is Gateaux- differentiahle in U. If 
the map: 

F^:W ^C{X,Y), F^ia)^d'^F 

is continuous at a* € W , then F is Frechet- differentiahle at a* and its Frechet derivative eval- 
uated in a* results 

do^.F^d^F 

Lemma 11 The map V defined by gip is C^{H^O,T); H^{0,T)). 

Proof. The continuity of V has been discussed in Lemma [6l Next we consider the differential 
map V : H^{0,T) C{H^{0,T), H^{0,T)). Due to Theorem \TU\ it is enough to prove that 
a -> d'^V is continuous w.r.t. the operator norm in £{H^ {0,T), H^{0,T)). The action of d'^V 
on u e i/i(0,T) is 

d^V{u)^q, (93) 
q = -u^-a(^uiO)e^'^? g^i;0)l^^ + ^Uq^ . (94) 

Then, the iJ^-bound 

(95) 

and the continuous dependence from a follows from Lemma [B] and proceeding as in Corollary [T] 
■ 

Proposition 12 For e D{Ha(Q)), the map a -> r{a) is of class C^{H^{0,T), H^il)). More- 
over, the regular part of the solution at time T, r{a) — V{a){T)QQ , considered as a function of 
a, is of class C\H^{0,T), f] H^{I)). 
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Proof. The continuity of T{a) follows directly from the continuity of the map a — >■ V{a) 
and the estimate ([77)1 in Lemma |51 Next consider the differential map F' : H^{0,T) — 
£{H^ (0,T), Hq{I)). Our aim is to prove that V is continuous in the operator norm. Due 
to the Theorem [TOl it is enough to prove that a — > d'^T is continuous. The action of d^F on 
u e H^{0,T) is: d^F(u) = F{d^Viu),T); for a, a e H^{0,T) the difference: d^T - rfgF is 
expressed by 

d^T - d^F - F(d^^(w) - d^V{u), T) 

Using ((77)) . one has 



|cr-4F|||= sup ||i^(d^i^M-d^^M,r)||^ < 

nei/^(0,T) ° 



< sup \\d'^V[u)-d'iV[u)\\^^ ^\\\d^V-d^V\ 

'ueH^{0,T) 



then, the continuity of a — > d^F follows from the continuity of a — >■ d'^V proved in Lemma [TTl 
Let introduce the map 

M(a)=F(a)-y(a)(r)go^(-,0) (96) 

representing the regular part of the state at time T. Combining the result of Lemma [TT] and the 
first part of this proof, one has: n{a) S C^{H^{0,T), Hq{I)). To achivc our result it remains 
to prove that: -dlfj,{a) belongs to C^{H^{0,T), L'^{I)). According to ([83l)-(l84l), this can be 
written as 

\e'^^'4>^^+F{w,T) - w{T)g^{-,0) + [V{a)] {0)Q{T) 



d^ 

dx^ 



f,ia) - -~f^[iF [w, T)-F {w, T)) - {w{T) - w{T)) G^] + 



where = [V{a)] (t) - [V{a)] (0), while the regularity of the terms: -^^e^^^^ cj)^ and Q{T) 
is discussed in Proposition [51 For a,a £ i?^(0,T), the difference (M(a) — f^{ct)) is 

dx^ 

-{[V{a)] (0)-[y(a)] (0))^Q(T). 

Using ([57)) . one has 

^i{a) - fiia) = i [F (w', T)~F (w', T)) + A (w(T) - w{T)) + 

-{[Via)] {0)-[Via)] m^Q{T) 

w' and w' respectively denoting the time derivatives of the functions V{a) and Via). The 
continuity of ^(a), then, follows from the continuity of F{-,T) (expressed by ([77|) and of the 
map V{-). The Gateaux derivative d^/i acts on u G H^{0,T) as 

^"^("^ ^ [^^(d«u;(^.),T) - [d'^wiu)] (r)C?o\.,0) + [d^Via)] (O)Q(T)] . 

Proceeding as before, the difference d^fj.{u) — d%ii{u) is 
d^M(") - rff m(") = » {F (dtd'^wiu) - dtdiw{u),T)) + A ([d^w;(u)] (T) - [df u;(u)] (T)) ^0^+ 

- ([d^u;(z.)] (0) - [d?z«(u)] (0)) ^Q(r) . 

Using dni) and recalling, from LemmafTTl that a d^V is continuous in the C (i?^(0, T), i?^(0, T)) 
operator norm, this relation leads to continuity of a — d^/i in the C (iJ^(0, T), L^(/))-norm. ■ 
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4.2 The linearized map 

Next we discuss the surjectivity of the hnearized map d'^T for a £ Hg{0,T). Since the trans- 
formations V and d'^V preserves Dirichlet boundary conditions, in this framework the result of 
Lemma [TT] raphrases as: 

Lemma 13 The map V defined by fHP is C\H^{0,T); H^{0,T)). 

This also implies that the state at the initial and finial times possesses only regular part. 
Thus, Proposition [T^ become: 

Proposition 14 For tpo e H'^r\H^{I), the map a r{a) is of class C^{H^{0,T), H'^r\H^{I)). 

It is whorthwhile to notice that, when ipQ is given by a linear superposition of eigenstates of 
odd kind (i.e. tjjk — sin |x, k even), the source term in ([?T|) . e**^i/'o(0), is nuU at each t, and the 
charge, V{a) in the above notation, is identically zero. In these conditions, the particle does not 
'feel' the interaction and the evolution is simply determined by the free propagator e**^? . This 
implies: T{a) — e'^'^'^' Vo- In order to avoid this situation, we assume the initial state tpo to be 
fixed in the subspace of even functions 



<y5 = ^ CfcV'fc ? (97) 
keD J 



This choice also implies that F takes values in W, as can be cheked by its explicit formula. Due 
to Proposition [Til one has: P £ C^{H^{0,T), W). 
The linearized map d^P is defined according to 

d^P(«) ^^Y.[[ li')^'''" ^"^"'^ A (98) 



keD 



- q{t) = -a{t)Uq{t) + u{t) e'*^Vo(0) + -U [V{a)] (t) (99) 

TT \ TT / 

Let ijjT S W be a target state for the linear transformation (|M)) - (|99p . We address the following 
question: does u £ Hq{0,T) exists such that d^r{u) = ^ 
Consider at first the equation 



^T = ^y^\ r p{s)e-'^^^^-'^ ds\ (100) 

w.r.t. the variable p. Denoting with the Fourier coefficients of ipT in the basis {'/'felfegi); the 
above equation is equivalent to 



Cfe = ^ r pis^e-'^"^^-'^ ds, keD (101) 
V"" Jo 

This is a trigonometric moment problem; the existence of solutions is established by K. Beauchard 
and C. Laurent in [5] for arbitrary values of T > according to the following result. 
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Theorem 15 Let T > and {ujk}i^f^f^^^, with No = NU {0}, be an increasing sequence of [0, +oo) 
such that ujQ — and 

hm (a;„+i - tj„) = +oo 
For every c G £^(No, C), there exist infinitely many w G ((0, T) , C) solving 

T 



Cfc= / v{s)e"^-'ds, nGNo. (102) 
Jo 

Proof. The proof rephrases the one given in Corohary 1 of [S], where the authors discuss the 
existence in i^((0,T),M) under the particular condition: c G ^^(No,C), cq G M. For sake of 
completeness we give a sketch of it. 

Let define: a;_„ = —ujn for n G Nq. According to the result of [13] (referred as Theorem 1 in 
[5]), the family {e*'^"*}^^^ is a Riesz basis of the subspace defined by the clousure in L^{0, T) 
of the Span This condition is equivalent to the existence of an isomorphism J (we 

refer to Proposition 20 in [8]): 



v{s)e^--^ ds 



J: J" — > ^{'L.C) 

V 

Given c G ^^(No, C), we consider c G £"^{1, C) such that: Cfc — Ck for fc G No- To any choice of c 
there corresponds an unique solution to (|102p defined by J~'^{c). ■ 

When the coincide with the frequencies of the standard basis {e*""*, uj = 
L^{0,T), equation (|10ip can be also interpreted as a constraint over the Fourier coefficients of 
the function p. In this case, it is possible to give a solution of the problem respecting Dirichlet 
boundary conditions. These remarks are summarized in the following proposition. 

Proposition 16 The following properties hold: 

1) For {ck} G £2(0) and T > 0, U01\) admits infinitely many solutions p G L^{0,T). 

2) Let T > 8tt. For any xpx G VV, it is possible to find at least one p G L[q{0,T) solving U0(^} . 

Proof. 1) Let Vt = 'l2keD '^ki'k, and consider the corresponding moment equation (llOip . For 
fc G I?, we set: fc = 2n — 1, a;„ = X2n-i with n G N, and loq = 0; with this notation, our problem 
writes as 

i-T 

p(s)e''""'* ds, n G No (103) 



P = A=P, c„,T|„gN = e'^'''^C2„-i and cq = Q (104) 

According to theorem [TSl (jl03p admits infinitely many solutions p G L^(0,T), which are deter- 
mined by fixing the extensions of the families e"^"'' and c„^t to n G Z. 

2) ForT — Btt, the functions {^~^^''^] ^eo form a subsystem of the standard basis je*'^"*, uj — 
in £^(0, T). Thus, the equation (|10ip partially determines the Fourier coefficients of the function 
p. A particular solution po is obtained by a superposition respecting Dirichlet conditions on the 
boundary 

^Po(i) = V e'^-^Cfc (e-^^'=* - e^^*-*) . (105) 



/TT 
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Since the regularity of e W imphes: { k'^Ck} G £2(0), it follows p e H^{0,T). For T > Stt, 
the target is attained by p such that: pI^^^q g^^j = Po is a solution of (jlOOl) in -ffo(0, Stt) and 
p(i > Stt) = 0. ■ 

In order to solve the inverse problem related to (P5 1) - (j99p . we proceed as follows: for a target 
state ipT S W, Proposition [TBI allows to determine (at least one) q g Hq{0, T) solving ([98]); then, 
this function is replaced in (|99p in the attempt of finding a suitable u E H^{0,T) such that 

- q{t) = -a{t)Uq{t) + u{t) f e**^^o(0) + -U [V(a)] {t)\ (106) 

TT \ IT J 

The solvability of this equation w.r.t. u is strictly connected with the specific choice of the initial 
state, V'o, and of the linearization point a. In particular, if the function (e^'^V'olO) + \V{a)\ (t)) 
is different from zero at each t, one can easly determine u as 

u{t) = (^-q{t) - ^a{t)Uq{t)^ [e^'^MO) + [V{a)] (t)^ (107) 

Consider a simplified framework where ipo is a Laplacian's eigenstate in W 

ipo=iPk k odd (108) 
and a ~ 0. In this setting, the equation (l99|) writes as 

- qit) ^ ^u{t) e-'^f'' (109) 

and the inverse problem 

A ( i-T \ 



is solved in iJQ(0,T) by the last point of Proposition [1^1 This result, and the regularity of the 
map F leads to the local steady state controllability of the system (j90p - (PT|) around the point 
a = 0. 

Theorem 17 (Local controllability) Let T > 8tt and assume il08]) . It exist an open neigh- 
bourhood V X P C Hq{0, T) X W of the point (0, ip^^) such that: T\y : V ^ P is surjective 

Proof. As remarked at the beginning of this section, for -00 e W, the map F is {Hq (0, T) , W) . 
Moreover, under our assumptions, c?qF is surjective for a = 0. Then, by the Local Surjectivity 
Theorem (e.g. [TTl) we know that there exists an open neighbourhood V x P C Hq{0,T) x W 
of (0, F(0)) such that the restriction r\y : ^ -J> P is surjective. By definition, F(0) = e~*^s^-(/'fe 
and P is a neighbourhood of i/'fc in W. ■ 
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